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Using density functional theory we have performed theoretical investigations of the electronic 
properties of a free-standing one-dimensional organometallic vanadium-benzene wire. This system 
represents the limiting case of multi-decker V n (C6H 6 ) n +i clusters which can be synthesized. We 
predict that the ground state of the wire is a 100% spin-polarized ferromagnet (half- metal). Its 
density of states is metallic at the Fermi energy for the minority electrons and shows a semicon- 
ductor gap for the majority electrons. We found that the half-metallic behavior is conserved up to 
12%, longitudinal elongation of the wire. However, under further stretching, the system exhibits a 
transition to a high- spin ferromagnetic state that is accompanied by an abrupt jump of the magnetic 
moment and a gain of exchange energy. 

PACS numbers: 71.20.-b, 75.50.-y, 75.75.+a, 85.75.-d 



During the last years, molecular magnets have been 
attracting enormous attention, because they are consid- 
ered as potential candidates for future applications in 
high-density information storage and quantum comput- 
ing. Among such novel systems, we would like to focus 
on vanadium-benzene multi-decker clusters (V n Bz n+ i, 
Bz=C6H6). These are molecular magnets appearing in 
a gas phase 0, 0. The synthesis of the V-Bz clusters 
was realized during a reaction of laser- vaporized metal 
atoms with benzene in He atmosphere |2, |3| . The mass 
spectrometric measurements jj, |^ and theoretical stud- 
ies [M El suggested that V n Bz n+ i clusters, approach- 
ing a length of up to n = 6 (18 A), should have a 
one-dimensional structure which is settled by the metal- 
ligand interactions confined to a single molecular axis. 

In this Letter, we present the study of the infinite one- 
dimensional (ID) organometallic vanadium-benzene (V- 
Bz) wire shown in Fig.l which could be considered as 
the limiting case of the V n Bz n+ i clusters with n — > oo. 
We have performed ab initio calculations of the elec- 
tronic structure of the system and predict that the wire 
under consideration is a new unusual example of the 
one- dimensional half- metallic ferromagnet. It has inte- 
ger magnetic moment of per unit cell and the spin- 
polarized band structure with finite density of states at 
the Fermi level for one spin channel and a semiconductor 
gap for the other spin channel (Fig. 2). Thus, the elec- 
tronic properties of this organometallic wire are similar 
to that of the known Heusler compound NiMnSb which 
has been predicted by de Groot et al. in 1983 to be a 
half- metal 0. Due to a spectacular electronic structure 
such wire or its finite fragment would act as a nearly 
perfect spin-filter with very high magnetoresistance ratio 
when it is placed between ferromagnetic electrodes in the 
two-terminal spin-polarized- STM like set-up 0. 

Our calculations are based on density functional the- 



ory and were performed within the linear combination 
of atomic orbitals formalism implemented in the highly 
accurate CRYSTAL code 10]. The Perdew- Wang [nj 
parametrization for the exchange and correlation energy 
was used. The description of the electron subsystem was 
done with the full electron basis set 6-311G** [13 for 
carbon and hydrogen atoms and with a frozen core ba- 
sis set SBKJC for vanadium. The convergent re- 
sults were achieved using the grid of 24 /c-points in the 
ID Brillouin zone. The wire under consideration (Fig.l) 
has a one-dimensional hexagonal lattice with P6/mmm 
(D\ h ) rod symmetry. From calculations of the total en- 
ergy of the system we have found that the optimal struc- 
ture corresponds to the following interatomic distances: 
dvBz=1.69 A, rcc = l-45 A and r C H = L09 A. These data 
are in a good agreement with other theoretical studies of 
V-Bz clusters [lJ,LL5|- We have also examined another, 
staggered- type, geometry of the wire when the consec- 
utive benzol rings were rotated by 30° with respect to 
each other. We have found that the staggered geometry 




FIG. 1: (Color online) The structure of the V(C 6 H 6 ) wire. 
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FIG. 2: (Color online) 
The spin-resolved band 
structure (left plot) and 
the DOS (right plot) of the 
V(CeH6) wire in the ferro- 
magnetic phase. The la- 
bels nearby the circles at 
the band structure plot in- 
dicate crystalline orbitals 
of the wire calculated for 
the r point. These func- 
tions are shown in the cen- 
ter. The corresponding or- 
bitals of the minority and 
majority spin bands have 
practically the same shape. 



is energetically less favorable, with a total energy being 
0.06 eV higher than the lowest-energy of the wire with of 
D\ h symmetry. Moreover, our calculations showed that 
the rotation of the benzene ring does not change the elec- 
tronic properties of the wire. 

In Fig. 2 we present the spin-polarized band structure 
of the V-Bz wire (left plot) together with the correspond- 
ing density of states (right plot). We accompanied these 
data by plotting the crystalline orbitals evaluated at the 
T-point. The symmetry of the bands corresponds to the 
Cq v point group relevant for the wire k- vectors. The 
notations of Ref.[l6| are used throughout the text. We 
see from Fig. 2 that for the minority spin electrons the A\ 
band and the double degenerate E2 band cross the Fermi 
level, while for the majority spin electrons the Fermi en- 
ergy lies in the gap. The direct gap has a width of 1.18 eV 
while the spin- flip gap is 0.12 eV at the T point. Such 
system would have ballistic conductance of 3 e 2 /h in one 
spin channel and zero transmission for the other spin. 
To verify the obtained result we have performed addi- 
tional calculations within the Vosko-Wilk-Nusair LDA- 
type 17] and hybrid B3LYP exchange-correlation 
potentials together with all electron [HI and frozen core 
LANL2DZ |20( basis sets. We have confirmed that the 
half-metallic behavior is robust against details of the cal- 
culations. The exact value of the gap width and relative 
positions of the bands are of course subject to the chosen 
functional. The LDA predicts a bit smaller value for the 
gap (around 1.0 eV) while the hybrid B3LYP functional 
noticeably pulls the bands apart the Fermi level. We have 
also examined the antiferromagnetic (AFM) state of the 
wire and found it to be energetically less preferable as 
compared with the ferromagnetic configuration. 

We show in Fig. 3 the total valence charge density 
(pj + pi) and spin density (p^ — pi) contour plots for the 
ferromagnetic V-Bz wire. The density plots are given for 
two planes: along the wire and perpendicular to it. The 
electronic charge density in the V-Bz bond critical point 



is 0.132 e/ajj. According to the Bader analysis |2l| of the 
topology of the charge density, the covalent type of bond- 
ing is predominant. The vanadium atom is charged neg- 
atively (— 0.22e) due to a charge transfer to benzene that 
agrees with the findings of Refs.|d l22j. Thus, the ionic 
type of bonding is also present. The magnetic density 
map (Fig. 3b) indicates a quite localized positive mag- 
netic moment at the V atom (+1.28/ie) and a small neg- 
ative magnetic moment (— 0.28/ie) redistributed over six 
carbon atoms. Our value of the magnetic moment at the 
V atom agrees well with the theoretical data of Ref.|(| 
obtained for finite V n Bz n+ i complexes where these val- 
ues varied from 1.15/i# up to 1.36/i# depending on the 
complex size and the vanadium position in the complex. 
The total magnetic moment of the wire unit cell is equal 
to 1.0 /jlb- The same integer value of the total magnetic 
moment was obtained in EPR measurements 113, 0, [2^ 
as well as in the calculations performed for the single V- 
Bz and multidecker V n Bz n+ i clusters with n < 6 ja.[2^. 

The analysis of the bands shown in Fig. 2 together with 
data presented in Fig. 3 leads to the following picture 
of the formation of the ground state. The fully occu- 
pied bands below —6 eV with respect to Ep arise from 
the to pure benzene orbitals which are not spin-splitted 
and thus play no role in the formation of the magnetic 
state. For energies above — 6 eV the wire bands are 
formed from the hybridized states of the vanadium and 
the benzene, i.e. from molecular orbitals of the V-Bz clus- 
ter. The strong crystalline field splits the vanadium 3d 
states to the singlet A\ state (d^ z i_ r i) and two dou- 
blets, of Ei (d xz , d yz ) and E2 (d xy ,d x 2_ y 2) symmetry. 
The doubly degenerate E\ band just above —6 eV is 
formed mainly from the HOMO and HOMO— 1 7r-type 
orbitals of benzene with an admixture of d xz , d yz vana- 
dium states. Due to the strong hybridization effect be- 
tween the states of the same symmetry, the vanadium 
d X z, d yz levels are pushed well above the Fermi energy 
where they are coupled with antibonding benzene states 
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FIG. 3: (Color) Total valence charge density (a) and spin density (b) for V(CeH6) wire. 



leading to the formation of the two different bands of 
Ei symmetry, marked by (f) and (g) in Fig. 2. The re- 
maining d 3z 2_ r 2 and d xyi d x 2_ y 2 vanadium states form 
two bands nearby the Fermi level, of A\ and E2 symme- 
try, respectively [they are marked by (b) and (c)]. These 
bands are spin-split. The vanadium atom itself has three 
3d electrons. It turns out that in the wire the vanadium 
s states are shifted above Ep and are responsible for the 
formation of the A\ band labelled as (e) in Fig. 2. Thus, 
two s electrons of vanadium move to the 3d shell and five 
electrons in total wish to occupy three levels {A\ and E2) 
which are available per spin. Finally, majority spin elec- 
trons fill completely two bands of A\ and E2 symmetry 
which therefore are placed below Ep . The remaining two 
electrons of each unit cell are redistributed among minor- 
ity spin A\ and E2 bands both of which are crossing the 
Fermi level. To conclude, the magnetic moment per unit 
cell is forced to be 1.0/i#. The semiconductor gap in the 
majority spin channel is formed between two different E2 
bands one of which originates mainly from d xyi d x 2_ y 2 
vanadium states while the upper one comes from LUMO 
and LUMO— 1 antibonding 7r*-states of benzene. 

Let us move again to Fig. 3 and consider the spin den- 
sity plot. We see from Fig. 2 that the E2 minority spin 
states nearby Ep are almost filled, therefore the mag- 
netic moment at the vanadium atom arises mainly from 
the A\ states [labelled as (b)]. As it is seen from Fig. 2 the 
corresponding wave function (b) is nearly a pure d 3z 2_ r 2 
vanadium orbital so that its symmetry is reflected in the 
angular distribution of the spin density around the V 
atom. We note further more, that the small negative 
magnetic moment on the C atoms results from the un- 
compensated E2 band (c) whose molecular orbital con- 
tains an admixture of the carbon p states (Fig. 2). 

We can ask further the following question: what is the 
origin of the exchange interaction between quite local- 
ized magnetic moments of V atoms? According to the 



structure of the wire, we can assume that different mech- 
anisms of the exchange are possible: direct exchange as 
well as the indirect, superexchange interaction which in 
the present case could be established through the 7r-type 
orbit als of benzene. To demonstrate that direct exchange 
between vanadium moments indeed takes place we com- 
pare the spin-projected local density of states (LDOS) at 
the V atom for the ferromagnetic (FM) and the antifer- 
romagnetic (AFM) state (see Fig^J). If there would be no 
mutual influence of localized moments, one could expect 
the LDOS for both phases to be nearly indentical. As it 
is clearly seen from Fig. 4, this is not the case for the V-Bz 
wire. The structure of the LDOS changes for the AFM 
phase. The AFM phase is not any more the half- metal 
and is characterized by the reduced magnetic moment at 
the V atom (0.62/i^). The direct exchange interaction 
between V moments favors the FM phase to be energet- 
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FIG. 4: (Color online) Local density of states at the V atom 
for FM and AFM configurations of the wire (in the AFM 
state the second V atom has inverted LDOS concerning spin 
direction compared to the first one). 
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FIG. 5: (Color online) Outer plot: Binding energy Evbz be- 
tween V and Bz as a function of the unit cell length. Black 
and red lines correspond to the low- and high-spin FM states, 
respectively. Blue line with open circles refers to the AFM 
state which for any lattice constants is less preferable as com- 
pared with FM state. Innner plot: Behavior of the magnetic 
moment per unit cell under elongation of the wire. 



ically preferable. However, we can not exclude another, 
indirect type of spin coupling, which is governed by the 
Kramers- Anderson mechanism [23, E3 of the superex- 
change interaction. For the considered system, because 
of the overlap of the HOMO (HOMO-1) 7r-type orbital of 
benzene with the d xz , d yz orbitals of vanadium, one of the 
p electrons from benzene can hop to one of the V atoms 
and the remaining unpaired p electron on Bz can enter 
into a direct exchange with the other V atom. Applying 
the semi-empirical Goodenough-Kanamori rules |22| re- 
garding to the sign of the effective exchange integral, we 
can assume that the indirect exchange interaction tends 
to aline vanadium moments antiparallel and thus it com- 
petes with the direct coupling. 

We have also studied the effect of the longitudinal 
stretching on the properties of the VBz wire. Our re- 
sults are summarized in Fig. 5 where we show binding 
energy between V and Bz as a function of the unit cell 
length (outer plot) together with the behavior of mag- 
netic moment (innner plot). We have found that a half- 
metallic state characterized by the integer magnetic mo- 
ment is conserved up to ~ 12%, elongation of the wire. 
At the critical unit cell length of about 3.75 A, the sys- 
tem exhibits transition to the high-spin state. Just before 
the critical point the partially occupied minority spin A\ 
band (b) nearby Ep is transformed to the very narrow 
peak pinned at the Fermi level. Under further stretching 
this band is pushed above Ep and a fractional charge of 
about 0.3e per unit cell from the minority spin band flows 
to the upper lying majority spin band of E\ symmetry 
(f) which now intersects the Fermi level close to the A 
point of the Brillouin zone. Thus, the transition is ac- 
companied by an abrupt jump in the magnetic moment 
(inserted in Fig. 5), from \\ib up to ~ 1.6/i# per unit 
cell, and high spin state becomes energetically favorable 



because of the gain in the exchange energy. 

To conclude, using the ah inito density functional 
theory based calculations we predict, at least theoret- 
ically, a new class of half- metal ferromagnets: infi- 
nite one-dimensional organometallic vanadium-benzene 
[V n (C6H6) n +i (n — > oo)] wires. We mention, that ac- 
cording to our calculations, the half-metallic properties 
have been also found for similar systems, where V was 
substituted by Mn or Co as well as for hybrid systems 
(Mi-Bz-M 2 -Bz) n (n -> oo), where Mi and M 2 stand for 
different (V, Mn or Co) atoms. 
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